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071 Atmospheric Oscillations. 
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TiitwcUidion, 

1. The chief question discussed in this paper (§§ 6 — 12) is that of the free 
oscillations of an atmosphere whose temperature varies v^ith the altitude ; and 
in particular the case of a uniform vertical temperature-gradient is studied in 
some detail. For consistency it is assumed that the expansions and con- 
tractions follow the adiabatic law. The problem is treated as a two- 
dimensional one, the space co-ordinates involved being horizontal and vertical ; 
and the more definite conclusions arrived at relate to the case where the 
(horizontal) wave-length is somewhat large in comparison with the height 
of the atmosphere. 

The results are most easily interpreted when the temperature -gradient does 
not fall much below that characteristic of a state of convective equilibrium. 
The normal modes of oscillation then fall into well-defined types. 

In the most important type, the motion of the air-particles is mainly 
horizontal, and independent of the altitude, and the waves may therefore 
be described as " longitudinal." The velocity of propagation of progressive 
waves is found to be equal to ^{gH), where H denotes what may be called 
the "virtual height" of the atmosphere, i.e. the height of a "homogeneous 
atmosphere " corresponding to the temperature of the loivest stratum. That 
the result should come out intermediate in value between the velocity of 
sound in the lowest stratum, viz., v^(7^H), and the zero velocity corresponding 
to the zero temperature which is postulated in the higher regions was to be 
anticipated ; but that it should be identical in form with that obtained on 
the hypothesis of an isothermal atmosphere whose expansions are subject to 
Boyle's law,* the effect of the upward decrease of temperature being exactly 
compensated by the greater elasticity implied in the adiabatic law, is 
somewhat remarkable. 

When the temperature-gradient falls distinctly below the "convective'' 
value, the character of the oscillation is less simple. The w^ave- velocity is 
somewhat increased, but must always remain below the value ^(7^H), which 
is the velocity of sound in the lowest stratum. 

2. A second type of oscillations depends on the degree of stability of the 
atmosphere. 

* Rayleigh, *Phil. Mag/ (4), 1890, vol. 29, p. 173 ; ^Scientific Papers,' vol. 3, p. 335. 
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The work required to bring unit mass of air from the density pi to the 
density p2 under the adiabatic condition j?//?'y=;^i/piy is 

Hence if we imagine two thin strata of equal mass, whose densities are pi, p2, 
and pressures pi, p2> to be interchanged, the work required to effect this will 
be, per unit mass, 

{7~l)pi IW J (7~l)i02LW J 7-1 W /02V- 

(2) 

If we avail ourselves of the notion of ** potential temperature/'* i.e, the 
temperature ^ which any particular portion of air would assume if brought 
adiabatically to some standard density D, we have 

p/py ^m/J)y~\ (3) 

where E is the constant of the formula 

P = -Rpe, (4) 

6 denoting the absolute temperature. Hence (2) becomes 

E(p2y-^-piy-^)(^l-~^2) /.x 

Hence if p2>pi, we must for stability have ^i>%; ie. the potential 
temperature must increase upwards. Now, if y denote depth below a 
standard level, we have, in equilibrium, 

dp/dy = gp ; (6) 

and combining this with (3) and (4), we find 

Id^ _1 f de {j-l)g' 



^dy evdy E 



(7) 



In convective equilibrium, where p/pv, and consequently d, is the same at all 
altitudes, we have 

dy'^ yR ' ^^ 

This equilibrium, though stable for some types of disturbance (§8), is in 
other respects neutral. For complete stability, d^jdy must be negative, and 
therefore 

f^krz^. (9) 

dy 7E 

When this condition is fulfilled, we have a series of possible modes of 

^ V. Bezold, ' Berl. Sitzb.,' 1888, vol. 46. 



1910.] On Atmospheric Oscillations. 553 

oscillation whose periods, depending as they do on the extent to which the 
temperature-gradient differs from the convective value (8), are comparatively 
long. Oscillations of this character, governed by local conditions, must 
undoubtedly occur in the atmosphere, and may conceivably account for some 
of the minor fluctuations of the barometer. 

There remains a third type of oscillations which, when the wave-length is 
moderately great, approximate to the character of waves propagated vertically 
in the atmosphere. These have been discussed in a previous paper by the 
author.'*' From a meteorological standpoint they can hardly be of importance. 

3. The theory of the " longitudinal '^ waves is of interest in relation to the 
large-scale oscillations of the earth's atmosphere as a whole. This subject 
was treated by Laplace,t and is of some importance in connection with the 
suggestion put forward by Lord KelvinJ as to the origin of the semi-diurnal 
variation of the barometer. Laplace's investigation was based on the 
hypotheses of a uniform equilibrium temperature and an isothermal law of 
expansion, and on the further assumption that the vertical motion of the 
ciir-particles may be neglected. § Since the circumstances are then practi- 
cally those of sound-waves propagated horizontally, his results naturally 
involve the " Newtonian" velocity of sound, ^(^H), where H is the height 
of the homogeneous atmosphere corresponding to the assumed uniform 
temperature do, viz., H = l^Oajg. 

The hypotheses referred to were, of course, adopted only for mathematical 
convenience. As representations of actual conditions they are very imperfect ; 
and there is, moreover, great uncertainty as to the most suitable value to be 
attributed to 6 a. It appeared to the writer that a firmer ground for quanti- 
tative conclusions would be gained if it were possible to calculate the 
wave-velocity (for long waves), even in the two-dimensional problem, on 
somewhat more natural suppositions as to the constitution of the atmosphere 
and the law of expansion. 

In the actual atmosphere the temperature, as a rule, diminishes upwards, 
although (as we have seen) it is necessary for stability that the gradient 
should nowhere exceed the convective value. The special hypothesis of 
a uniform gradient, which is here adopted as a basis of calculation, is itself an 
artificial one ; but in spite of the fact that it implies an upper limit to the 
atmosphere, it may claim to give, on the whole, a better representation of the 

* 'Lond. Math. Soc. Proc' (2), 1908, vol. 7, p. 122. 

t ' M^canique Celeste,' Livre 4, Chap. 5. See also Eayleigh, loc, cit, 

% *Eoy. Soc. Edin. Proc.,' 1882, vol. 11 ; * Math, and Phys. Papers,' vol. 3, p. 341. 

§ Some such assumption is necessary to make the problem determinate, in the absence 
of a prescribed condition to be fulfilled, or approximated to, \yl the upper regions of 
the atmosphere. 

2 Q 2 
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true conditions than the isothermal view, on which, indeed, the earth's 
atmosphere is merely a local concentration of a medium diffused through 
space. 

As regards the law of expansion, since permanent inequalities of tempera- 
ture are postulated in the equilibrium condition, it is proper to ignore the 
transfer of heat between adjacent portions of the air during the oscillations. 
In any case, theory shows that the effect of conduction on such long waves as. 
we have here in view may safely be neglected.* 

The main conclusion of Laplace was that the free and forced oscillations of 
an atmosphere covering a globe, whether this be at rest or in uniform 
rotation, are identical with those of a liquid ocean of uniform depth H ; but 
in view of the nature of his premises, and of the uncertainty as to the 
temperature to be adopted in estimating the value of H, considerable doubt 
has been felt as to how far this analogy can be relied upon for qnantitative 
resnlts. The present investigation tends, I think, to show that inferences of 
this kind will not be very far from the truth, provided the temperature 
adopted be the mean temperature of the lower strata of the earth's atmo- 
sphere, so far as this can be ascertained. The formal adaptation of the theory 
of longitudinal waves to the case of an atmosphere of relatively small depth 
covering a globe would follow the same course as in Laplace's investigation. 

4. As regards the semi-diurnal variation of the barometer, the passage of 
Kelvin's paper already referred to runs as follows : — 

" The cause of the semi-diurnal variation of barometric pressure cannot be 
the gravitational tide-generating influence of the sun, because, if it were, 
there would be a much larger lunar influence of the same kind, while in 
reality the lunar barometric tide is insensible or nearly so. It seems, 
therefore, certain that the semi-diurnal variation of the barometer is due to 
temperature. Now, the diurnal term, in the harmonic analysis of the: 
variation of temperature, is undoubtedly much larger in all, or nearly all, 
places than the semi-diurnal. It is then very remarkable that the semi- 
diurnal term of the barometric effect o/the variation of temperature should be 
greater, and so much greater as it is, than the diurnal. The ' explanation 
probably is to be found by considering the oscillations of the atmosphere, as 
a whole, in the light of the very formulae which Laplace gave in his 
* M^canique Celeste ' for the ocean, and which he showed to be also applicable 
to the atmosphere. When thermal influence is substituted for gravitational, 
in the tide-generating force reckoned for, and when the modes of oscillation 

* This follows from the equations (due substantially to Kirchlioff and Rayleigh) given 
in the author's ^ Hydrodynamics,' 3rd edit., § 343. Radiation has a different tendency in 
this respect. 
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corresponding respectively to the diurnal and semi-diurnal terms of the 

thermal influence are investigated, it will probably be found that the period 

of free oscillation of the former agrees much less nearly with 24 hours than 

does that of the latter with 12 hours ; and that, therefore, with comparatively 

small magnitudes of the tide-generating force, the resulting tide is greater in 

the semi-diurnal term than in the diurnal." 

The first question which here arises, viz., whether as a matter of fact the 

earth's atmosphere has a mode of oscillation of the requisite type, with a 

period of about 12 mean solar hours, can at the present time be examined 

more closely than was possible at the date (1882) of the above extract. The 

free oscillations of an ocean of water of uniform depth covering a globe of the 

size of the earth, rotating with the same angular velocity, have been very 

fully investigated by Hough* in the course of his classical work on tidal 

theory. He finds, in particular, that in the case of the most important free 

oscillation having the same general character as a semi-diurnal tide wave 

{i.e. its most salient spherical harmonic constituent is the sectorial liarmonic 

of the second order), the depth h for which the period is exactly 12 sidereal 

hours is given by 

gh/4:(o^a^ = 010049, 

where a is the earth's radius, and « its velocity of rotation. This is evaluated 
at 29,182 feet. It is to be remarked, however, that throughout the calculation 
the mutual attraction of the disturbed fluid has been taken into account, 
whereas in the aerial ocean this influence must be quite insensible. If the 
disturbance were accurately of the type of a spherical harmonic of the second 
order the requisite modification would consist merely in multiplying the 
previous result by the factor 

^ = 1 — 3 X 0-18093 = 0-89144, 

where the decimal fraction in the second member is the ratio of the density 
of the water to the mean density of the globe, as adopted in Hough's 
computation. This would make 

ghjAw^a^ = 0'08958. 

As the result of a more direct calculation, using Hough's algorithm, together 
with such of his numerical results as are applicable, I find 

gh/4:(o^a^ = 0'08986, 
the last figure being somewhat doubtful. If we put g = 32*200, 
ft) = 2 7r/86164, 6^ = 20,902,000, this gives 

h = 25,930 feet. 

^ 'Phil. Trans.,' A, 1897, vol. 191, p. 139, See pp. 164, 179. 
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The substitution of a mean solar for a sidereal half-daj as the period 
involves a further slight diminution, which can be estimated pretty closely 
from another of Hough's results. He finds that for gh/4:0D^a^ =0*1, the speed 
(a) of the free oscillation in question is given by (t/co = 1*9968. Com- 
paring this with the former result, we infer that for a period of 12 mean 
solar hours (o-/co = 1*9945) we must have ghjAco^a^ = 0*09965, about. 
Assuming that when mutual attraction is ignored this figure is to be reduced 
in the same ratio as the former one, we have, finally, 

^A/4a)V = 0*08911, 

or, with the previous numerical data, 

A =25,710 feet. 

It must be remembered, of course, that these numerical results can claim 
no greater accuracy than the theory on which they rest, in which, in 
particular, the ellipticity of the earth, which is of the order 1/300, is 
neglected. 

On the other hand, the value of H for air at 0^ C. is about 26,200 feet, 
with an increase of about 96 feet for qyqtj degree above this temperature. 
The mean temperature of the air near the earth's surface is usually estimated 
at 15° C. This would make H == 27,640 feet ; but a somewhat lower value 
for the mean temperature of the lower strata, away from the immediate 
influence of the ground, would perhaps be more appropriate. 

Without pressing too far conclusions based on the hypothesis of an atmo- 
sphere uniform over the earth, and approximately in convective equilibrium, 
we may, I think, at least assert the existence of a free oscillation of the 
earth's atmosphere, of " semi-diurnal " type, with a period not very different 
from, but probably somewhat less than, 12 mean solar hours. 

At the same time, the reason for rejecting the explanation of the semi- 
diurnal barometric oscillation as due to a gravitational solar tide seems to call 
for a little further examination. The amplitude of this variation at places on 
the equator is given by Kelvin as 0*032 inch. The amplitude given by the 
'' equilibrium " theory of the tides is about 0*00047 inch.* Some numerical 
results given by Hough in illustration of the kinetic theory of oceanic tides 
indicate that in order that this amplitude should be increased by dynamical 
action some seventy-fold, the free period must suffer from the imposed period 
of 12 solar hours by not more than 2 or 3 minutes. Since the difference 
between the lunar and solar semi-diurnal periods amounts to 26 minutes, it 

■^ The numerical values given on p. 520 of the author's * Hydrodynamics ' relate to the 
lunar tide, and are, moreover, by an oversight, stated as *' amplitudes," instead of as 



" ranges." 
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is quite conceivable that the solar influence might in this way be rendered 
much more effective than the lunar. The real difficulty, so far as this point 
is concerned, is the a prm^i improbability of so very close an agreement 
between the two periods. The most decisive evidence, however, appears to 
be furnished by the phase of the observed semi-diurnaP equality, which is 
accelerated instead of retarded (as it would be by tidal friction) relatively 
to the sun's transit.* 

5. The concluding part of the paper (§§ 13, 14) is an attempt to examine 
more closely than has hitherto been done the theory of waves on a surface 
of discontinuity in the atmosphere. That such waves may play a part in 
meteorological phenomena has been pointed out independently by Helm- 
holtzf and Lord Kelvin,t but both writers have confined themselves to 
analogies drawn from the case of superposed homogeneous liquids. It is to 
be observed that even on this view the disturbance extends, upwards and 
downwards from the plane of discontinuity, through a space which is an 
appreciable fraction of the wave-length ; hence, apart altogether from the 
influence of compressibility, the conditions* of the question will be modified 
when the wave-length is such that the ordinary variation of density within 
this space becomes sensible. It seemed worth while to investigate the 
matter; but it must be acknowledged that when there are no currents, the 
discontinuity being one of temperature and density only, the analogy proves 
to be adequate, under such conditions as are likely to occur in the 
atmosphere, for a considerable range of wave-lengths. For very long waves 
it would break down, the disturbance ceasing to be even approximately 
concentrated in the neighbourhood of the plane of discontinuity. The 
discontinuity then becomes, in fact, an unimportant incident in the general 
upward diminution of density. 

When there is a discontinuity of velocity, the upper fluid being in steady 
horizontal motion relative to the lower, the question, when compressibility 
is taken into account, is more difficult, and I have not been able to arrive 
at any very simple results. There can be no doubt, however, that the 
aforesaid analogy is sufficient in this case also for wave-lengths less than a 
certain limit. In particular, the dynamical instability pointed out by 
Kelvin§ will remain. 

^ *Brit. Ass. Eep./ 1908, p. 606. The forced tides due to diurnal and semi-diurnal 
waves of temperature have been studied by Margules, 'Wien. Sitzb.,' 1890, vol. 99, 
p. 204. 

t * Berl. Sitzb.,' 1889 ; ' Wiss. Abh.,' vol. 3, p. 309. 

X ' Brit. Assoc. Rep.,' 1876 ; ' Math, and Phys. Papers,' vol. 4, p. 457. 

§ * Math, and Phys. Papers,' vol. 4, p. 76. 
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Theory of Long AUnospheric Waves. 

6. We consider an atmosphere arranged in horizontal layers of uniform 
density. The motions contemplated are restricted to two dimensions, x, y, of 
which X is horizontal and y vertical, the positive direction of y being down- 
wards. The equilibrium values of the pressure, density, and temperature 
are denoted by ^o^ po, -^o ; these are functions of y only, and are subject to 

the hydrostatic condition 

dp^ijdy ~ ^/>o, (10) 

as well as to the general relation 

po = E/)o^o. (11) 

The equations of small motion are, in the usual notation, 

cu dp dv dp , ,-ta\ 

where B/Bt = d/dt + ud/dx-i-vd/dy. (14) 

The expansions being supposed subject to the adiabatic law, w^e have also 

where c^ = ypo/po = T^^o^ 0-^) 

i.e. G is the velocity of sound corresponding to the equilibrium temperature at 
the point considered. It is accordingly in general a function of y. If we put 

j9=j^9o + C7, p = po+S, (17) 

and continue to neglect small terms of the second order, we have 

du dm dv djjj , ^ /i o\ 

9f+,^o^_ /|^+8£\ (19) 

ot ay \dx oyj 

Also, from (15), (13), and (10), 

dm . J)p (du , dv\ .^n^ 

5-+^M- = j^ = -7i^o(3-^+9-^j- (20) 

Hence, eliminating B and m, we find* 



3% ., 8 /9i6 . dv\ , 3 



c'~ ~+^ +g 



'V 



dt^ dxxd dyl dx' 

3% _ „ _3 (Su dv\ ,_^. (du dv\ dv 



y (21) 



y 



^ It may be noticed, parenthetically, that in the case of an isothermal atmosphere 
where c is constant, these equations are satisfied by 

u ~ e-(y-^) 9i//c^ f(ct - s), V = 0. 
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dv d' 



If we now write 

dx dy dx dy 

we deduce from (21), by differentiation, 



r, 



(22) 






(23) 



85 



(24) 



where y^ = d^/dx^-^d^Jdy^. (25) 

The latter equation shows that an irrotational motion is not possible 
unless 

(26) 



rf.c^ / 1N d0o (7~1)^ 



which we have seen to be the case of convective equilibrium. We note also 



that 






yPox> 



(27) 



-(7-1).} g. (28) 



hj (20). 

Eliminating f between (23) and (24), we obtain 

dt^ ^ df ^\ dy ^^^hydf ^ I dy 

If we assume that x and t occur only through a factor e^i^^-^^^)^ the 
equations (21) take the forms 

aht-i-ighv = —ikc^^, ^ 

9y 



— iglcu + a^v = — c^-?r^— 7^%, 



whence 






~\ 



y 



dy 



(30) 



^ 



Prom these, or from (28), we have 






'^'-^'-{^j-^^-M'^\ 



X 



= 0. 



(31) 



7. So far, the vertical distribution of temperature is arbitrary. In the 
case of temperature diminishing upwards with a uniform gradient, to which 
we now proceed, there is an upper limit to the atmosphere. If we take the 
origin of y at this level, we have 

^0 = ^y, (32) 
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where ^ is the gradient in question. It easily follows from (10) and (11) 
that 



)o cc ^», 



Pqoc;^ 



n+l 



where 

-owlSO 
jLxexxwe 



n = g/Ji/3 — 1. 
c2 = -yE/32/ = jgy/iw+l). 



(33) 
(34) 
(35) 



(cf*-9%2)t^ = ^ 



n+l 



^y^^ + (n+l)x-^^kyx 






\ {2'gJ + (^ + l)%l-'^2/% 



y' 



(36) 



and y||+(^+2)|^+|!i±i^%^I=±zl4_^a;fc^ = 0. (37) 

The meaning of the factor (ny— ti— l)/7, which appears in one of these 
terms, is to be noticed ; viz. we have 



7 ^ 



(38) 



where ^i is the temperature-gradient in a state of convective equilibrium, as 
given by (8). 

To solve (37) we put % = e^^^, (39) 



and obtain 



where 



^^+(^^+2 + 2%) ^ +2fltM = 0, 



(40) 



(41) 



This is integrable by series, the solution which is finite for ^ = being 

a . a+l 






a 



(2%) + 



(2%)2— ... ; 



(42) 



(43) 



(44) 



l.ti + 2 """'"^ ' 1.2. 1^1+2. 1?^ + 3 
or, in the notation of Dr. E. W. Barnes,* 

The remaining solution of (40) is of the form 

where <^ stands for the series in (42). This is not admissible in the present 

* See, for example, * Camb. Trans,,' vol. 20, p. 253, where references to other papers, 
are given. 

If we had assumed x == ^"^^ 4>^ in place of (39), we should have found 

<^ = jFj (^ + 2 — a ; ^ + 2 ; 2 %). 
The comparison verifies a well-known identity ; see Barnes, loc, cit 
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question, since it becomes infinite as y"^"^ for infinitesimal values of y, 
whereas the condition to be satisfied at the upper boundary is J)p/Dt = 0, or 
y«+i^ = ; see equations (27), (33). 
The formulae (36) now become 



71 -f 1 L dy 



Icycjy 



e^y, 



y (45) 



dy • ^- • V ^ ■ \ gjc^ 

the factor e^i*^^'^^^^ being omitted here, as elsewhere, for brevity. 

The condition that 'y = at the lower boundary, where y = h, say, taken 
in conjunction with (41), determines the values of a and a, the wave-length 
(27r/Jc) being supposed given.* 

8. In the case of oscillations about convective equilibrium we have 

;3 = /3i, 71 = 1/(7-^1). (46) 

It follows from (24) that d^^/dt^ = ; hence either f = 0, i.e. the motion 
is irrotational, or the period is infinitely long. 

The conditions to which the steady rotational motions thus indicated 
are subject follow most directly from (21). These equations are now 
equivalent to 



dx 






= 0, 



3 r 9/3^ . 8 



> 



< c 



dv\ 



+ ^ )+gv > = 0, 



(47) 



by (26). Hence 



dy L \3^ dy/ 

(T,_l)3,g + |) + , = const. 



(48) 



The choice of two functions, u, v, to satisfy this equation, together with the 
two boundary conditions, can be made in an infinite variety of ways. 

The remaining types of disturbance are periodic in character. The 
formulae (42) and (45) apply, with 



2 a 



gh 



+ ^?' + 2, 



(49) 



in place of (41). Since (42) makes 

(50) 

* The case of apparent failure, where 0-2 = gk^ does not arise. This would require, 
by (36), 

yg|+(^ + l-%)X = 0, or x = %-'*-i6^y, 
which violates the condition at the upper boundary. 
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the condition that v = for y = h may be written 

(% + l)2fiFi(a; n+1; -27ch)-(l-'QkhiFi(cc;n+2; -2kh) = 0. 
gk \ gkj 

(51) 

A complete discussion of the equations (49) and (51) is out of the 
question, but the limiting form to which the results tend as the wave-length 
increases is easily ascertained. In the first place, it appears that when kh is 
small we have 



(J 



2 



hh 



approximately, since this ensures, by (49), that cikh is also small. If H 
denote the virtual height of the atmosphere, as defined above (§ 1), we have 

H = h~^ [V% = h/in-^l), (53) 

Jo 

The limiting value of the wave- velocity V is accordingly given by 

V^ = 0-7P = ^H. (54) 

The bearing of this result has been discussed in the introduction. 
The formulae (45), (50), and (39) now lead to 

V oc k (y—h), J 

the factor ^^(0^^+^'^) being understood. These values fulfil, as they ought, the 
irrotational condition 

1=1= *"■ (=« 

Since the ratio of the amplitude of v to that of u is of the order hh, the 
motion is mainly horizontal, and the present type of waves may accordingly 
be characterised as '' longitudinal." 

The remaining solutions of (49) and (51), when hh is small, involve finite 
as distinguished from infinitely small values of ahh. As will be seen 
presently (§ 11), they approximate to the character of waves propagated 
vertically in the atmosphere. 

9. In the general case, where 71 is not restricted to the precise value 
lf{ry — l), the relation between a and a is as in (41). When hh is small we 
have still a longitudinal wave for which a^/gk is of the order M, subject to 
a certain condition. The equation (51) leads again to the result expressed 
by (52) or (54), and substituting in (41) we find that the implied assumption 
that ukh is also small will be justified provided /3i//3 — 1 be small, i.e. 
provided the temperature -gradient falls only a little short of the convective 
value ySi. 
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The limiting form of (42)^ when no assumption is made as to the order of 
magnitude of aM, is 

_. 2ahy (2aJcyf ' .^^. 

^ l.n + 2 l,2.n + 2.%^3 '"' ^ ^ 

or in Dr. Barnes' notation, 

^ = qFi (%'^2; —2aky\ (58.) 

whilst (50) becomes 

.2/7:^ + (ti4-l)<]b = (ti + l)oFi(^ + l; '—2uky). (59) 

oy 

It appears from (41), without making as yet any special assumption as to 
the smallness of ^i//3— 1, that when ahh is finite, whilst kh is small, the 
ratio cr^/gk will be very small or very great. 
In the former case we have 

$ = (f-')-t (^^) 

ultimately, and the condition (51) becomes 

(^+l)(|i-l)oFi(^^+l; -2^m)~2^MoFi(»+2; -2aM) = 0. (61) 

Since ©Fi (« + 1 ; - «) = H (n) ^"i^J^ (2z^% (62) 

in the notation of Bessel's functions, this may be written 

I «J.+i (m) = [j-l) J. (o>% (63) 

provided m^ = 8akh. (64) 

If ci) be a root of (63), the corresponding frequency of oscillation is given by 

jk'~'\w r~^' ^^^^ 

and the wave- velocity by 

V^ = (|-l).,M£:y)£l, (66) 

H being defined as before by (53). This result again is accurate as a limiting 
form for increasing wave-length. 

10. The equation (63) might be discussed, when n or 2n is integral, with 
the help of the tables of BesseFs functions, but it may be sufficient to 
consider the case where the ratio Oi— i8)//3 is small. It may be noticed 
that the formula embraces all the modes of the present class, the longitudinal 
waves already discussed corresponding to the case of c» infinitesimal. The 
roots of (63) which relate to the remaining modes are now given by 

J»+i (ft)) = 0, (67 
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approximately ; and in particular the first of these slightly exceeds the first 
finite root of (67). 

In convective equilibrium we have n = 2*55 if 7 = 1-40. The first finite 
root of JjpC®) = is G)i = 7, very nearly. Hence for oscillations about a 
state of very nearly neutral equilibrium we have 

V = (|j-l)*xO-53v^(^H). 

In the case of w = 3, which makes {/3i—^)/^ = ^, a first approximation, 

given by (67), is coi = 7'586, and a second is easily found to be c»i = 7*624. 

This leads to 

V = 0198v/(^H), 

which is about one-fifth the velocity of the longitudinal type of waves. 
As to the character of these slow rotational* modes, we find from (24) 






r=-^(§-l).r^.^^.%. (68) 



or, by (65X ^ = -^ • ^ • %• (69) 

Having regard to the kinematical meaning of the functions ^, x> ^s defined 
by (22), we see that the rotational quality in the relative motion of a fluid 
element predominates over the dilatational. We learn also from (45) that 
when 2/ = the amplitude of v is to that of tt in the ratio cr'^/gk, which is 
small. Since v vanishes at the lower boundary, we infer that the vertical 
component of the velocity is in general relatively small. The distribution 
of horizontal velocity depends ultimately on the function 

y 

which varies as {(oy^l1i^)~'^Sn{^y^P^^)> 

if CO be the relevant root of (63), or less accurately of (67). In the case of 
the first root, after the small one, this expression changes sign once, and 
once only, as y increases from to h. For n = 3, the change of sign occurs 
for (oy^lU = 6-379, or yjh = 070. 

The general character of the types of disturbance at present under 
consideration is most easily apprehended in the case of a " standing " 
oscillation. If on the preceding expressions we superpose others which 

•* The rotational character is, of course, present also in the longitudinal waves, unless 
^= jSj exactly, though to relatively slight extent. 
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differ only in the sign of <r, and reject the imaginary parts, we find, on 
discarding all but the most important terms, that 

30 



uoz 



2/ ^ + (^+1)0 I sin yfccc cos o-jf, 



I 



(70) 



is accordingly satisfied by 



'y X -—hy^ cos hx cos at. 
The differential equation of the lines of (oscillatory) motion, viz. 

vdx—%tdy = 0, 

y^'^'^f^^inhx = const., 

^i(w+l)J^^J(ft,yl/^l)siny^^ = C. 

If we put C = we get the lines y = 0, y = hy hx = stt^ but the former 
of these is only an approximation. The annexed figure indicates, without 
any attempt at minute accuracy, the general arrangement of the lines in the 
case of the lowest finite root of (67). 



or 



(71) 

(72) 
(73) 




In the modes corresponding to the higher roots there are horizontal nodal 
planes {v = 0), in addition to the lower boundary. 

Ee turning for a moment to the more important "longitudinal'' type of 
motion first considered, we note that the formulae (52), (54), cease to be 
accurate, even as limiting forms, when the ratio /Si//8 differs appreciably 
from unity. The formulae (65) and (66) will, however, still apply, cd being 
the lowest root of (63). As a numerical example, take the case where the 
temperature-gradient has only one-half the convective value, so that 

(/3i-/3)/^ = 1, ^ = 6. 

I find that the lowest root of 



IS (O 



: 4*96, approximately, whence 

V = l-07y(^H). 

The result must, of course, in any case be less than ^^(yg'H.), or 
1*18 ^(^H). The change of wave-velocity is accompanied by a change in 
the character of the oscillation, the variation of horizontal velocity with 
altitude now becoming sensible. 

The preceding formulae might also be used to estimate the rapidity of 
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falling away from the state of unstable equilibrium which prevails when 
/3 !> /3i, the value of a^ given by (65) being then negative. 

11. The modes for which a^jgh is large are easily accounted for. We 
have from (41) 

T = ^. (74) 

gk n^l 

and from (45) 3/^+(^ + l)^ = 0, (75) 

these being approximations which gain indefinitely in accuracy with increase 
of wave-length. On the present supposition that oikh is finite, notwith- 
standing the smallness of M, (75) reduces to 

J.(a)) = 0, (76) 

provided aP = SaJch, as in (64). If w be a root of this equation, the corre- 
sponding frequency is given by 



'yet' 



2 ___ yo) g ^^j^ 



4(71 + 1) * h' 

These modes are in the limit identical with the waves of vertical 
displacement discussed in a paper already cited in § 2. The formulae (45) 
show, in fact, that the ratio of the amplitude of i^ to that of v is for the 
most part of the order gJc/a'i If we put h = IjU^ti^ = iygTi^/{ri + 1), 
the equation (76) takes the form J^ ((rn) = 0, which is identical with 
equation (88) of the paper referred to. 

12. It may be worth while, for the sake of the contrast, to give the theory 
of the oscillations of a heterogeneous but i7icompressiUe fluid, whose 
equilibrium density has a similar distribution. 

Wehavenow Po^ = §' P«ar^-a^+^^' ('^> 

^ + A»=0, (79) 

at ay 

^+^=0. (80) 

OX cy 

If we put :p — po-i~m> p — /?o + S, (81) 

as before, we have 

?=-'«'^». (83) 

ot ay 



From (80) we have 

u 



_9t v = ^ (84) 

■by' dx' ^^ 
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i|r being the stream-function. Substituting in (82), (83), and eliminating m 
and S, we find^ 

|v^+l^( 3^+^11) =0. (85) 

If X and t occur only in the form e^(o-*+^'^^), we have 

p j^^+iplof+S!f^) = 0. (86) 

Also 'y = ^/^i|r, (87) 

and ^==^^+,^ko^i^^ i^^^ (88) 

If we now assume that po <^ y'^i (89) 

we have .^ J^^ + ^ |l^ + iVM^^hy] hf = 0, (90) 

or writing -s^r z= e^y<p, (91) 

2,^+(^ + 2%)|^ + ^(l+g)#= 0. (92) 

The solution which is finite for 3/ = is 



a /n7 \ . a . a + 1 



or </> = 1^1 (« ; '«' ; — 2%), (94) 

if 2« = 7^(l+2|V ^ (95) 

The second solution becomes infinite as y^^^ for 3/ = 0, and is therefore 
excluded, in virtue of (88), by the condition that Dp/Dt must vanish at the 
upper boundary. Since, by (87), i|r = for 3/ = A, we have 

iFi{a;n] -2]ch) = 0. (96) 

This determines a, and the value of a^ follows from (95). 

It is obvious that, when kh is small, (96) is not satisfied by finite values 

of a. If a be large, but so that akh remains finite, the equation (96) tends to 

the form 

ori(^^; -2aM) = 0, (97) 

or J,-i(ei)) = 0, (98) 

provided ©^ = Sakh (99) 

If ft) be a root of (98), we have 

^' = ^ = i| . kh, (100) 

gk 2 a €0^ \ / 

* Cf. Love, « Lond. Math. Soc. Proc. ' (1), 1891, vol. 22, p. 307. 
VOL. LXXXIV. — ^A. 2 E 
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and therefore, for the wave-velocity, 

V^ = ™ = ^''^''+^-^ . ^H. (101) 

Thus, if, for the sake of comparison with § 7, we put n = 2*5, we have 

> 

o>i = 4*491, whence 

V = r32v/(^H). (102) 

That the frequency should be increased by the incompressibility was to be 
expected ; that the effect is so considerable is due to the great modification 
which is caused in the character of the fundamental modes. 

The modes corresponding to the higher roots of (98) have horizontal nodal 
planes, and the frequencies form, by (101), a descending series,* as in the case 
of (65). 

Waves at a Surface of Discontinuity, 

13. When we proceed to examine the case of waves propagated along a 
horizontal plane where the equilibrium temperature is discontinuous, it may 
be sufficient to suppose the temperature uniform throughout each of the 
regions, above and below this plane, to which the influence of the waves 
extends. The plane in question is taken as the plane 3/ = 0, and the 
dependent variables relating to the upper region will be distinguished by 
accents. 

The formulae of §6 will therefore apply to the lower region, with the 
simplification that c is a constant; so that (31) becomes 

This is satisfied by ■)(^ = Ge^y, (104) 

provided \2+T|x+ X^'-F+fclk!^") = 0. (105) 

We are seeking for a type of motion analogous to that of waves on the 
interface of two liquids of different densities, in which case the values of X are 
Hh^. We assume, provisionally, that in our case also the roots of (105) are 
real and of opposite signs ; moreover, since the disturbance is to vanish for 
y z=z Qo , the negative sign is to be taken. 

For the upper region we shall have 

X = OVV (106) 

with a similar determination of V ; but the positive root is now the appropriate 
one. 

■^ Cf, Eayleigh, ' Lond. Math. Soc. Proc' (1), 1883, vol 14, p. 170 ; ^ Scientific Papers,' 
vol. 2, p. 200. 



1910.] On Atmospheric Oscillations. 569 

If 7] denote the ordinate of the surface of separation, as affected by the 

waves, we have 

dyjjdt = V (107) 

for 2/ = ; and the pressure in either fluid at the point {x, t]) is to be 
found by putting y = in the corresponding value of the expression 

i?o + tF + ^%. (108) 

Differentiating with respect to t, we see that Dp/Bt must be continuous at 
the interface * This involves, by (27), the continuity of %, so that the 
constants C, C in (104) and (106) must be equal. Again, it follows from 
(107) that V must be continuous, whence, by (30), 

(T^G^X--gkh^-\'yga^ = ah^^X' —gJ^c'^+fyga\ (109) 

This, together with the two equations of the type (105), determines the 
values of X, X^ and a. 

To obtain a solution, we denote the two equal members in (109) by fi ; 
thus 

X = 4,+-??-^, V = -^+^-^-'^. (110) 

Substituting in (105), we have 

yU.2 + ^(2yt;V-70-2) + (o-*-^2p)(o.2_p^)g2 = Q, (111) 

with a similar equation in which c is replaced by c'. Writing these two 
equations in the form 

ya2 + P/i+Q = 0, /i2+P>+Q' = 0, (112) 

and eliminating fi, we have 

(P-P')(PQ'-P'Q)+(Q-Q')' = 0. (113) 

Now 

P-P' = 2^F(c2-c'2;, •-) 

Q_Q' = {c^-e'^) {a^-g^]<?) {a^-W (c^ + c'^)}, i (114) 

PQ'-P'Q = ^(C2_C'2) (o-*-^F) {^o.4_^^2p (c2+ c'2)+ 2^*C2C'2}.^ 

Hence 

(<7*-/F){0-2-7c2(c2+c'2)}2 +2^2p{^o.4_^^a/^2(g2 + (,'2) + 2^*cV2} = 0. 

(115) 
This is of the fourth degree in a^, but one root only is relevant to the 
present question. The common root of (112) is 

^ Q-Q' ^^ ^a2_7,2(c2+c'2)' (lib) 

* This might almost have been assumed at once ; but it is to be observed that it would 
not give the correct condition to be satisfied at the common boundary of two currents. 

2 E 2 
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I jL 1 v/ V • JL ^x I 



whence 



X' _ ^^/<^ cr^ — M (r?'^ — e^) 



> 



J 



(117) 



&2 — -|(c^2— -C^X 

ft) = gjlm^^ 



If we now write 

the equation (115) becomes 

a?2(a?-2)2 + a)2{(27--.l)a.^^-4(7-.]l)a?-46V^*'} 

whilst - = -A ~~ «, 







X' _x—2Wlo/ 



It is to be noticed that 

52 



h x{x—2) 



(£f. 



(118) 
(119) 

(120) 
(121) 






r.2 



<f _ p— p 



,2 



a* c ^ + <? p + p 

where p, p', are the equilibrium densities at the plane y = 0, on the two sides. 
For sufficiently small wave-lengths, m and a? are very small, and the root of 
(120) with which we are concerned is x = mb^/a^, approximately; whence 



0-2 __ P P-gJc^ X = — ^, X^ = /v, 



(123) 



as in the case of superposed incompressible fluids.* 

To examine the matter further, the simplest procedure is to tabulate the 
function 

.,.2 x^(x'-'2y C124> 

4¥/a^~t4:{j-l)x-{2y'-l)x^ ^ ^ 

for a series of suitable values of x. The only case of real interest is where 
the discontinuity of temperature is very slight, so that h^/a^ is a small 
fraction. The following table gives a few results calculated on the supposition 
that P/a2 = yi^, with 7 = 1-40. The abrupt step in temperature then 
amounts to -^ of the mean of the temperatures (absolute) above and below. 



m. 



wK 



{JO. 



10A/(ip/ft>). 



-Xjh. 



VjJc, 



WaTe- 

length 

(metres). 



Period 
(seconds). 



10~« 
10-3 
10-4 

10"3 



-9960 X 10-8 
0-9615x10-6 

-^143 X 10-4 
2-004 xlO-3 



-9980 X 10-4 
-9806 X 10-3 
-8462 X 10"2 
-44^6 X 10-1 



1 -0005 
1-010 
1-088 
1-495 



0-998 


0-998 


0-981 


-980 


0-860 


0-841 


0-470 


0-426 



7*0 

69-2 

697-0 

8160 -0 



21-2 

65-9 

180-0 

301-0 



* Stokes, ' Camb. Trans./ 1847, vol. 8, p. 441 ; * Math, and Phys. Papers/ vol. 1, 
p. 197. 
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The fourth column gives the ratio of the frequency to that of waves of the 
same length on the surface of separation of two homogeneous liquids with 
the same discontinuity of density, as given by (123), viz., the ratio is 

on the present suppositions. The seventh column is calculated from 
^TTJk = 27ra?wlgy taking a = 332 metres per second ; and the last from 
27r/o- = 27raco/g^x, It is seen that, with increasing wave-length, the wave- 
velocity tends more and more to exceed the value estimated on the assumption 
of the homogeneity and incompressibility of the two fluids. At the same 
time, the disturbance tends to become, relatively as well as absolutely, less 
and less concentrated in the neighbourhood of the plane of discontinuity. 

14. In this question, again, it is of some interest to compare the case of 
waves on the common boundary of two liquids, each of v/hich, though incom- 
pressible, has a similar gradation of density. We therefore write, in (86), 

El = q, a. constant. (125) 

If we assume that 

i/r =Ce^y, (126) 

we derive \^ + qX+ (^^--l]k^ = 0. (127) 

These formulae may be supposed to relate to the lower region ; for the upper 
region we write q\ X\ C for q, \ C, respectively. 

The continuity of v involves, by (87) , that of yjr, so that C = C. Also, in 
virtue of the continuity of Dp/Dt, we have from (88) 

p (a'X +gk^) = p' ((T^V +^F), (128) 

where p, p\ are the densities just below and just above the plane y = 0. 

If the two fluids had been portions of the same gas at different 
temperatures we should have had 

q = yg/c^ q^ = jg/c'^ (129) 

and therefore p/p = q/q^. (1^0) 

Now from (127) we have 

q{a^\4-gk^) = a^{k^--\^). (131) 

Hence, if we adopt the relation (130) for the sake of the comparison, we 
must have \^ = X'^, or, taking account of the signs, X = —W This 
leads to 

k k q_-\-q o" 

and (^y-^_l^^£f-(ii::£;y = 0. (133) 
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The positive root of this quadratic in a^/gk is to he taken, since it is the 
only one which gives the proper signs to X, \\ it being assumed that />>p' 
and therefore q > q\ For infinitesimal values of q, q\ we reproduce the 
relations (123). 

In order to make the variations of density follow exactly the same law as 
in the atmospheric problem of § 13 we must give to q, q\ the values (129). 
In the notation of (118), (119), we have then 

^^^r^m^x-m%^la^ = 0. (134) 



x^ 



The following table, like the former one, refers to the case of 
V^ja^ =1/100. In order that the comparison may be for the same series of 
wave-lengths, those values of m are chosen which were obtained in the previous 
numerical work. The significance of the column headed 10^(^/(») is the 
same as on p. 570. The comparison shows the usual effect of a constraint in 
increasing the frequency. 



w. 


10^{xl<v), 


K'lh, 


Wave-length 
(metres). 


Period 
(seconds). 


-9980 X 10-4 
-9806 X 10-3 
0-8452x10-2 
-4476 x 10-1 


1 -003 
1-035 
1-324 

2-534 


0-993 
0-934 
0-570 
0-156 


7-0 

69-2 

597-0 

3160 -0 


21-2 

64-3 

148-0 

178-0 



